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Table 1 

Matrix elements of local and conserved vector currents. The quantity My depends on the coefficient C 
in eq. (9). Some common factors are suppressed, but the third column gives the bottom line. 



matrix element 



local 



conserved 



should be 



{Q(p)\Vo\Q(p)) 




ZyG 


1 


{Q{pm\Q{p)) 


Zy^e-^^yi/My 


Zya])ilM2 


Vil'mq 


{0\Vo\Q{p)Q{-p)) 











{omQ{p)Q{-p)) 


Zy^e-^^"" 

I- 


ZyaCXl + \CE sinh(2Mia)]e-*fi" 


1 



level 



X = 2M2a 



C(l + moa) 



C 



moa[2 + rrioa) 

Finally, G '^^ — Gi and A'^' — A are suppressed by 
a power of a, so I neglect them here. 

If c_B = (so Mb ^ M'j), M'j/Mb 
as nioa co; hence the Wilson action un- 
derestimates the chromomagnetic contribution to 
Am- If C = (instead of di, eq. (7)), X ^ 2C 
as nioa co; hence the usual local current 
overestimates the current contribution to Am- In 
both cases the relative error in Am is 0(0*^). In 
contrast, maladjustments in the heavy-quark cb 
and di produce no error in cfjoo ■ 



On the other hand, if 



CB(moa,gQ) 



di{moa, gg) can be computed reliably in tadpole- 
improved perturbation theory, (our version of) 
the improvement program reduces Am — A^m ^° 
0((/o')- Nonperturbative adjustments of and 
My could reduce the difference further. 



5. CONCLUSIONS 

In our analysis of the action [2,3,6], we found 
that there were only two large effects associated 
with the large mass limit of actions for Wilson 
fermions. One is the distinction between M2 and 
Ml , and the other is the field normalization factor 
gMia/2 [1-3^6]. Eqs. (6), (8) and sect. 3 show how 
the mass dependence of the unobservable field 
normalization filters down to currents and other 
multi-quark operators. 

With mass-dependent improvement the re- 
maining lattice artifacts are 



SO ~ dn\moa) X (ap)dim 0„-dim 



where p is a dynamical scale. Experience, both 
perturbative [4] and nonperturbative [9], indi- 

(13) 

cates that the residual mass dependence interpo- 
lates smoothly from small-mass to static results, 
with a "knee" at M2 between and 5. In partic- 
ular, the c„ and C„ do not grow uncontrolledly. 
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and spin-blind Coulomb term generates the only 
dynamics. The corrections to the flavor symme- 
try are correct only if M2 = niQ (so tune the bare 
mass mo a), and the corrections to the spin sym- 
metry are correct only if Mb = rnq (so tune the 
coefficient of the clover term [7] cb)- 

From the heavy-quark effective theory [8], one 
knows to construct currents, etc., not with the 
non-relativistic held, but with 

^>i{x) = [l-fD/{2mQ)]^'nn{x)- (5) 
Combining eq. (4) and (5) 

*i(a;) = exp(|Mia) [1 + adif-D] ij{x), (6) 
where 

1 C(l + nioa) 1 



*1 — £,NR 



2mQa mo a (2 -|- mo a) 2M2a 



(7) 



The last equality assume M2 = mq to obtain the 
correct amount of heavy-quark flavor symmetry 
breaking. With ^i{x) operators take the simple 
form 



jI^(x) = *((a;)r*i'(a; 



3. IMPROVED OPERATORS 



(8) 



Consider the construction in eq. (2). There is 
a redundancy because, for example, one can use 
either local or point-split operators. Consider a 
modiflcation of the local current 

= Zy.^[j, + aC{yD 7, + j.fD}]^. (9) 

The role of the dimension-4 term will emerge 
shortly. Table 1 contains matrix elements of 
and the conserved gauge current . {V^ differs 
from the Noether current if the action has a a ■ F 
term.) All three nonzero matrix elements are cor- 
rect for the local current, if Zyh = Zyh = e^'-" 
and if C is tuned so that My = M2 ■ 

The "form-factor" matrix elements of the gauge 
current would require Zya — Zya — 1, but 
the "decay-constant" matrix element would re- 
quire a different Zya for Mi ^ 0. One can ob- 
tain all three matrix elements correctly only by 
taking Zya ^ Zya and adding to a term 
aC'ipDitp, as in Vl" . But the new term spoils 



the attraction of the gauge current, namely that 
current conservation determines its normalization 
exactly. Hence, the gauge current seems to be 
even less practicable than the local current. 

To tune My = M2 at tree level one must choose 
C = di. Indeed, with ZyL = e^^" , eq. (9) is the 
same as the construction of eqs. (5)-(8). One 
flnds that rfi « (1 — (^Q)moa/4: for small mass, 
when the only 0(a) correction is in Zyh. On the 
other hand, one flnds that di = (2mQa)~^ in the 
inflnite mass limit. Without the improvement, 
the operator is incorrect at 0(l/mQ). 



4. ERRORS ANALYSIS FOR 



'M 



Improved Wilson fermions obey heavy-quark 
symmetry, even when the quantities cb and di 
are maladjusted. Thus, one can exploit the ma- 
chinery of heavy-quark effective theory to esti- 
mate the associated errors. We shall illustrate the 
method with the decay constant (fiM = \/mM /m , 
where M = (P, V) denotes pseudoscalar and vec- 
tor heavy-light mesons. Similar results hold for 
semi-leptonic form factors. 

According to heavy-quark effective theory 



1 



A 



M 



niM 



(10) 



where cfjoo is the matrix element, evaluated at in- 
flnite heavy-quark mass. The l/niM correction 
has three distinct contribution, 



Am — Gi + QhMG2 + |/imA, 



(11) 



from the heavy-quark kinetic energy, chromomag- 
netic interaction, and l/mq expansion of the 
weak current (eqs. (5) and (8)) [8]. The helic- 
ity factors are hp = 1 and hy = —1/3. Each 
contribution to Am can be related to a matrix 
element in the inflnite mass theory. [Gi > and 
A = niM - mq.) 

With eqs. (3) and (9) the slope becomes 



(12) 



where the second line assumes mq = M2 ■ At tree 
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We discuss lattice artifacts for matrix elements of hadrons containing one or more heavy quark. In particular, 
we analyze interrelations between lattice artifacts and the l/mg expansion. The implications for calculations of 
heavy-light decay constants and of semi-leptonic form factors are discussed. 



1. INTRODUCTION 

During the past several years, we have dis- 
cussed how to make sense of Monte Carlo compu- 
tations of Wilson fermions, even when the fermion 
mass becomes large [1-4]. (By "Wilson fermions" 
I mean a four-component field and Wilson's so- 
lution [5] of the doubling problem.) Because the 
numerical data never satisfy mqa <^ 1 and of- 
ten have AqcD Z^TiQ ~ |-, we choose not to as- 
sume either limit when calculating couplings of 
the improved action, or renormalization of com- 
posite operators. If warranted, these limits can 
be recovered afterwards. 

This report concentrates on external operators, 
whose matrix elements are needed to study elec- 
troweak properties of heavy-light hadrons. For 
comprehensive exposition with lattice terminol- 
ogy, see ref. [6]. To complement that article, 
this one often uses the language of heavy-quark 
effective theory instead. Nevertheless, eqs. (7) 
and (13) are valid for all mqa and mg/AqcD [6]. 

The completely renormalized (i.e. cutoff-free) 
action is 



5(mQ/AqcD) = ^CnimQa,gl)Sn, 



fl) 



where the sum runs over all interactions Sn with 
the desired field content and symmetries. Simi- 
larly, the completely renormalized operator is 

= Zo {{moa},gl) ^ C„{{moa},gl) 0„ , (2) 



where the 0„ are all operators with the cor- 
rect quantum numbers. In practice one truncates 



the sums and computes couplings c„ , normaliza- 
tion factors Zo , and coefficients C„ only approxi- 
mately. We match on-shell amplitudes computed 
via tadpole-improved perturbation theory in gg. 
At every fixed order in gg, however, we match 
to all orders in mo a. (Higher-order corrections 
Cn"' are mass dependent too, as clearly stated in 
ref. [3] and explicitly demonstrated in ref. [4].) 

Here we truncate eqs. (1) and (2) at tree level 
(/ = 0), dimS'n < 5, and dimO„ < 1 + dimC 
Hence, our construction, like any other, yields an 
action S = S + SS and operators O = O + SO 
that are imperfect. To make sense of a numerical 
calculation with S and O, one must understand 
the errors stemming from SS and SO. 

2. HEAVY-QUARK SYMMETRY 

Treating the gauge field semi-classically, the 
Hamiltonian (logarithm of the transfer matrix) 
can be written [2,6] 

i2 



H = ^ 



NR 



Mi+joAo 



2M2 



iS-B 

2M, 



B 



NR 



(3) 



The non-relativistic field ^nr(2^) is related to 
tp{x), which appears in the action, by a Foldy- 
Wouthuysen-Tani transformation and a certain 
normalization factor [2,6]. Neglecting terms of 
order (pa)^ and higher, 

'^nR(x) = exp(^Mia)[l + a^nRfD]ij(x). (4) 

The "masses" Mi , M2 , and Mb depend the cou- 
plings C, mo a and cb, defined in ref. [2,6]. 

The Hamiltonian obeys heavy-quark symme- 
tries for infinite mass. In that limit the fiavor- 



